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We study a slow classical system [particle] coupled to a fast quantum system with discrete energy 
spectrum. We adiabatically exclude the quantum system and construct an autonomous dynamics 
for the classical particle in successive orders of the small ratio e of the characteristic times in order 
to uncover new physical phenomena. It is known that in the order e° the particle gets an additional 
[Born-Oppenheimer] potential, while in the order it feels an effective magnetic field related to the 
Berry phase. In the order the motion of the clasical particle can be reduced to a free [geodesic] 
motion on a curved Riemannian manifold, with the metric generated by the excluded quantum 
■ system. This motion has a number of unusual features, e.g., it combines subspaces of different 

^— ■») ' (Riemannian and pseudo-Riemannian) signature for the metric tensor. In the order the motion 

, of the classical particle is still described by a Lagrangian, but the latter linearly depends on the 

particle's acceleration. This implies the existence of a spin tensor [non-orbital angular momentum] 
for the particle. This spin tensor is related to the momentum via an analogue of the zitterbewegung 
effect. The Hamiltonian structure of the system is non-trivial and is defined via non-linear Poisson 
brackets. The linear dependence of the effective classical Lagrangian on higher-order derivatives is 
' seen as well in the higher orders e" 
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I. INTRODUCTION 



' A recurrent theme in modern physics is to derive autonomous equations of motion for an open system, i.e., a system 
^ ' that interacts with its environment [1-3]. Depending on the type of environment, there are different conditions under 
O^. which this procedure is possible. A group of methods, which goes under the name of system-bath interaction, amounts 
to isolating a relatively small system in contact to an equilibrium environment (thermal bath) [1, 2]. The prerequisite 
, of applying the model reduction in this case is that the reaction of the system on its environment is in a sense weak 
J> ■ [1, 2]. One of the main consequences of this approach is the Langevin equation, which supplements the Newton 
\^ \ equation of motion for the small system by a (random) conservative and and non-coservative (i.e., non-Lagrangian) , 
0^ . velocity-dependent friction force [1, 2]. 

There is another set-up that allows deriving autonomous equations for an open system. Here the essential condition 
is that the target system is much slower than its environment [3-9]. One of the oldest results in this direction is 
l/^ the Darwin Lagrangian [3, 4]: when the characteristic speed v of charges is slow as compared to the speed c of the 
, electro-magnetic field, the latter can be adiabatically excluded, producing at the order {v/cY the Darwin Lagrangian 
for the charges, which has important applications in plasma physics and astrophysics; see [4] for a recent review. 

In this paper we shall study the quantum-classical (also called mean-field or hybrid) dynamics, which describes 
coupled quantum and classical systems [5, 6]. This is the most used set-up for coupling quantum and classical 
variables [10], and has numerous applications, e.g., in chemical physics [11] and in (semi)quantum gravity [12]. 
Assuming that the classical system is slow — a condition that is normally met in practice — we exclude the fast 
' quantum system and study to which extent the ensuing dynamics of the slow classical system can be described by an 
autonomous Lagrangian-generated equations for the classical coordinates. For example, in the above Darwin problem, 
the autonomous dynamics exists up to higher orders in (w/c), but already at the order {v/cY the dynamics of the 
charges is not Lagrangian due to the friction (Abraham-Lorentz) force related to radiation damping, i.e., to emission 
of electromagnetic field by the charges [3]. (In some special cases this friction force is suppressed and the Darwin 
Lagrangian can be written including the order (w/c)^ [3]). 

For the quantum-classical dynamics it is well known that at the zero order of e — where e is the small ratio of the 
characteristic times for the quantum over the classical system, respectively — the influence of the quantum system 
on the classical one can be described by the Born-Oppenheimer potential energy term [5-7, 10-12]. It was shown 
by Berry and Robbins that in the first order of e one gets an effective magnetic field, which manifests itself as the 
velocity-dependent term in the classical Lagrangian [5]. Goldhaber has recently shown that in the second order 
one gets in the Lagrangian an additional kinetic energy term, i.e., a quadratic form in slow velocities [6]. A very 
similar result on the order was obtained earlier by Weigert and Littlejohn for two coupled (fast and slow) quantum 
systems [7]. 

What happens in the next orders? In particular, how far we can continue the expansion over e, still keeping the 
classical system Lagrangian? Most importantly, are there new physical effects essentially related to post-adiabatic 
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corrections? 



Here we answer these questions. It appears that at every order over e one can derive Lagrange equations for the 
dynamics of the classical system. However, there is an important difference between the orders e and and all 
successive orders. At the order the classical dynamics is Lagrangian, but the Lagrangian starts to depend on the 
higher-order time-derivatives of the classical coordinates: While the classical Lagrangians normally depend on the 
coordinates and their first-order time-derivatives (velocities), at the order we get a Lagrangian that is linear over 
the classical accelerations. 

This fact is of conceptual relevance. The classical physics is essentially based on the Newton's second law that 
equates acceleration to the force, which depends only on coordinates and velocities. As a consequence, the trajectory 
of the classical motion is fixed via initial coordinates and initial velocities. In its turn, the Newton's second law is 
generated by a Lagrangian, which depends on coordinates and velocities. A Lagrangian depending on higher-order 
derivatives enlarges the amount of the initial data needed to fix the classical trajectory and produces equations of 
motion that go beyond the Newton's law. Such Lagrangians were phenomenologically introduced at various places 
and for various purposes (see, e.g., Ref. [13] for a sample of references), but our result seems to be the first example 
where a higher-derivative Lagrangian emerges for an open classical system due to time-scale separation. We should 
like to stress that the fast system being quantum is not important for obtaining the above result. What is important 
is that the fast system is integrable, i.e., it admits a full and globally well-defined set of action-angle variables. 

Dependence on higher-order derivatives in the Lagrangian implies a number of essential changes in the kinematics 
of the classical system: the momentum of the classical system depends on the acceleration, while the full angular 
momentum tensor is a sum of the usual orbital part and a term that can be interpreted as the spin of the classical 
system. In the simplest non-trivial case this spin is proportional to the velocity square of the classical particle. 
We show that this implies the existence of the zitterbewegung effect, where the momentum of the classical particle 
(system) is governed by the projected time-derivative of the spin. So far the zitterbewegung effect was known only 
in the physics of relativistic Dirac electron; see [14] for a review, while we show the same effect appears in a purely 
non-relativistic slowly evolved classical system due to its coupling to a fast quantum system. It appears now that this 
effect is a part of the physics generated by higher-order post-adiabatic corrections. Similar dependence on higher-order 
derivatives is expected at higher orders e" with n > 4, though in the present paper we restict ourselves with deriving 
the effective classical Lagrangian up to the order e^. 

While these results concern higher-order (three and more) post-adiabatic corrections, we found an interesting effect 
already in the second-order post-adiabatic correction, which was formally known since Ref. [6, 7]. It appears that 
the slow classical motion within this order can be reduced to a free motion on a Riemannian space with a signature- 
indefinite metric tensor. This implies a possibility of interchanging between time-like and space-like coordinates. 
Recall in this context that within non-relativistic classical mechanics the geodesic motion on a curved surface proceeds 
according to a positively-defined metric tensor, while the geodesic motion in the general theory of relativity has a 
metric tensor with signature (1, —1, —1, —1) [3]. In both cases the signature is fixed. 

This paper is organized as follows. In section II we introduce the quantum-classical dynamics. The next section, 
section HI outlines the adiabatic perturbation theory, which differs from the standard text-book presentations by a 
careful accounting of the higher-order terms. In section IV we review the derivation of the classical Lagrangian in the 
orders e and e^. In particular, we reproduce in a systematic way the results obtained by Berry and Robbins [5] and 
Goldhaber [6]. At the order this classical Lagrangian corresponds to a classical particle moving along the geodesies 
of a curved manifold. We calculate the curvature for the simplest non-trivial case and work out its implications for 
the stability of the effective classical motion at the order e^. Here we also point out at an unusual scenario related to 
the metric of the manifold changing its signature [i.e., chaning from a Riemannian to a pseudo- Riemannian manifold]. 
Section V derives the classical Lagrangian at the order e^, and shows that the classical Lagrangian in this case depends 
linearly on the third-derivative of the classical coordinates. In this section we explore kinematical consequences of this 
result and explore its Hamiltonian description. In section VI we deduce the classical Lagrangian at the order and 
show that it also depends linearly on higher-order derivatives of the classical coordinates. The last section presents 
our conclusions and summarizes the present work. Several technical question are discussed in appendices. 



II. QUANTUM-CLASSICAL DYNAMICS 



Consider a K-degiee of freedom classical system with coordinates q= {qi, . . . ,qK) and with Lagrangian 




(1) 



where M is the mass, and V{q) = V{qi, . . . , qx) is the potential energy. 
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Now this classical system (or particle) couples to a quantum system with Hamiltonian operator H[q{t)], which 
parametrically depends on the classical coordinates. The quantum system evolves in time according to the Schrodinger 
equation (for simplicity we put h= 1) 

idt\^) = H[q{tmit)), (2) 

where is the wave- function, and where 3t = 

The classical part of the dynamics is written as [10-12] (see [10] for a derivation) 

M^+d^V+{'^{t)\d^H[q{tm{t)) = 0, ii = l,...,K, (3) 



dt2 



where we defined ^ : 



Eq. (3) is the Newton equation of motion, where besides the classical force —d^V, the classical particle experiences 
an average force —{^{t)\d^H[q(t)\\^{t)) exerted by the quantum systems. In this sense the classical coordinates play 
a role of a mean-field [10]. The main purpose of the present paper is to understand to which extent this force can be 
generated by a Lagrangian which depends on the classical coor(linatc;s qa and their time-derivatives. 

It should be clear from (2, 3) that the total average energy is conserved in time: 

^{t)\H[q{tm{t))^ =0. (5) 
We note that the quantum-classical equations of motion (2, 3) can be derived from a Lagrangian 




c = - - (mm + f Ell ( 



dt 



V{q), (6) 



where as a set of independently varying parameters one should take j^) and q (or alternatively and q)"^. It is seen 
that £ is simply a sum of the corresponding quantum and classical Lagrangians, which points out at the naturalness 
of quantum-classical equations of motion (2, 3). 

Let us briefly comment on derivations of the quantum-classical dynamics from a full quantum-quantum dynamics. 
Such a derivation was carried out in literature several times; sec, e.g., [10, 15-17]. Moreover, many derivations of the 
(semi) classical mechanics from the quantum mechanics can be adopted for deriving quantum-classical dynamics; see 
in this context [18, 19] in addition to the above references. The main assumption involved in all these derivations is 
that fluctuations of classical coordinate(s) are small [10, 15-19]. For deriving the quantum-classical dynamics it is 
not necessary that the classical motion as such is slow ^. Note that the derivations of the quantum-classical dynamics 
need not neglect fluctuations of all pertinent variables, i.e., it need not impose the full quantum trajectories. It will 
suffice that the to-be classical sector of the dynamics is approximated via suitable Gaussian density matrices [18]. 
Then, the parameters of this matrices satisfy the equations of motion for some efii'ective classical systems [18] In 
Appendix A we briefly remind the main argument involved in the derivation of quantum-classical dynamics. 



^ Note that 9^ = 9qj,(t) acts only on the coordinates, but not on the velocities, e.g., duqa = 0. In particular, 9^ commutes with the total 
time-derivative 4-. 

^ As usual, when varying (6) we put aside the total time-derivatives, e.g., ^((5'J'|*). 

^ It is not excluded that there are situations, where both the classical limit and adiabatic limit — where the classical motion is slow — are 

taken simultaneously. Our consideration docs not apply to such situations, because when looking for post-adiabatic corrections one 
should simultaneously account for post-classical corrections, which is something we do not do. 
* A general remarks is in order here. The quantum-classical dynamics is just an approximation which holds under suitable conditions. 
There are, however, controversial aspects related to this dynamics, which emerged when people wanted to get a non-perturbative 
generalization of the mean-field quantum-classical dynamics; see, e.g., [20] for examples. These generalizations are supposed to be closed 
and self-consistent theories, where one part of variables is quantum and another is classical. Such theories (if they exist) would somehow 
get the same fundamental status as their limiting eases, i.e., as quantum and classical mechanics. Numerous attempts to formulate 
such fundamental quantum-classical theories met with severe difficulties [20]. Those diSiculties do not seem to be insurmountable, as 
witnessed by a recent proposal by Hall and Regginato [21]. 
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A. Classical representation for the quantum system 

Below we are going to concentrate on the adiabatic limit of the quantum-classical system, where the classical 
system is slow and the quantum system is fast, and derive an autonomous cqiiations of motion for the elassical part. 
A natural question is that why specifically we need the fast system to be quantum, and would it be possible to obtain 
the same result assuming that also the fast system is classical. (Then we would not need additional conditions for 
the applicability of the quantum-classical dynamics, and we could start from the outset with an overally classical 
dynamics) . 

The answer to this question is that in principle only one feature of the fast quantum dynamics is needed, that 
is its integrability [in the sense of [22]]. To support this answer we may note that all the results of the quantum 
adiabatics (including the quantum adiabatic theorem and the Berry phase) have their natural counterparts in the 
classical adiabatic theory for integrable system (i.e., the Berry phase becomes the Hannay angle in the classical theory) 
[5]. 

There is a more direct way of seeing the analogy between the quantum dynamics and the classical integrable 
dynamics [30]. The Schrodinger equation (2) can be mapped to a classical dynamics if one introduces a classical 
Hamiltonian 

n = imi'f) = E ^nmKTm, (7) 
Knm = (7n|^l7m), = (7m I*) (8) 

where {|7n)} is some time-independent orthonormal base in the considered Hilbert space. Now introducing the angle 
{4>n) and action (J„) variables via the phase and modulus of r„ as r„ = ^JT^e^^'^ ^ one can show that the Schrodinger 
equation for r„ reduces to the classical Hamilton equation 

d<f)n din 

This is a classical integrable dynamics, because it has globally well-defined action-angle variables with trivial Poisson 
brackets [22]. Eq. (9) makes clear that the adiabatic theory for the quantum Schrodinger equation can be alternatively 
developed from the viewpoint of classical integrable systems [30] . 

For definitness, in the present paper we shall confine ourselves with the quantum fast system. 



III. TIME-SCALE SEPARATION AND ADIABATIC PERTURBATION THEORY 



We shall now assume that there is a time-scale separation: the quantum system evolves much faster than the classical 
particle. To make this assumption more precise and to investigate its consequences, let us recall that the adiabatic 
energy levels {Ek{q)}t=i and the corresponding eigen-vectors {\k{q))}'^^i are defined via the eigen-resolution of the 
Hamiltonian H[q] at fixed values of g = {qi, . . . ,qK)- 

H[q]\k{q))=Ek{q)\k{q)), {k{q)\l{q)) = 6ki, k = l,...,d, (10) 

where d is the total number of energy levels. We shall assume that the adiabatic energy levels are not degenerate. 

Then qualitative sufficient condition for the time-scale separation is that the characteristic time of the classical motion 
is much larger than where A is the minimal adiabatic energy gap: A = mmk^i{\Ek — Ei\). ^ 
Note that the adiabatic representation (10) has a gauge freedom: 

|fc(g))^e^«.(«)|fc(5)), (11) 

where ak{q) is an arbitrary single- values function of q — {qi, . . . ,qN)- Hence all physical oservables have to be 
gauge-invariant . 

To reflect mathematically the fact of time-scale separation we shall write the dependence of the quantum Hamilto- 
nian on the classical coordinates as 

H[qi{et),q2{et),...], (12) 



This condition is sufficient, but not necessary for the validity of the time-scale separation and the consequent adiabatic approach, e.g., 
the latter can still hold if certain level-crossings are allowed. We shall not consider this more general situation in the present paper. 
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where e is a small dimensionless parameter 

e«:l. (13) 

The time-scale separation, i.e., condition (13), can be generated, e.g., by a large mass M of the classical particle. 
Then the classical particle moves slowly — provided that its initial velocity is small — and e ^ l/\fM . This scenario of 
time-scale separation is normally met in chemical physics (heavy classical nuclei versus light quantum electrons) [11] 
and semi-quantum gravity [12]. 

In the Schrodinger equation (2) we shall assume that the initial state |^'(0)) is equal to an eigenstate: 

|*(0)) = |n(0)). (14) 

Within the adiabatic approach the choice (14) does not imply any serious loss of generality; see Footnote 6. 

Our program is now to solve the Schrodinger equation (2) imdcr the adiabatic assumption (12, 13), and determine, 
via this solution, the structure of the averaged force in (3). To this end, we shall need the adiabatic perturbation 
theory, which was developed in [23, 24], and which is explained in detail in Appendix B. Now we shall recall some 
basic facts from this theory. As in any theory that is based on time-scale separation, we should start with dividing 
the sought solution into fast and slow components: 

\ii{t)) = e-'J'o^"^"''") |^„(ei)) (15) 
= e-i/o<i"^="(")|V;„(s)), (16) 

where in (16) we introduced the slow time s = et, and where [■^„) satisfies [see (2)] 

ie\Ms)) = [H{s) - £„(s)]|V„(s)). (17) 

Here dots denote differentiation over the slow time 

s = et, A=—, (18) 

and the lower index n in (15, 16, 17) refers to the initial state (14). Depending on the context we shall write H[q{s)] 

as H{s), etc. 

Eqs. (15, 16) extend to the adiabatic situation the usual formula for a stationary state of a time-independent 
quantum Hamiltonian. This analogy also explains why \tpn) in (15, 16) depends only on the slow coordinate. We 
see that the dynamical phase e~i -l^o duE„(u) jg ^j^g jg^g^ component of the wave- vector, since due to e -C 1 it stronlgy 
oscillates at slow times ^. 

Within the adiabatic perturbation theory, the solution of (17) can be sought for via expanding over e (see [23, 24] 
and Appendix B) 

|V'«(s)) =e^"'^" |0„(,s)), (19) 

|<^„) = \n{s)} + e|ni(s)) + e>2(s)) + e'|n3(s)) + • • • (20) 

The zero order term |(/)„) = |n(s)) in this expansion is the statement of the adiabatic theorem. In (19), ("(")!"(")> 
is the Berry phase factor; it was separated out for ensuring the proper gauge-covariance [24]; see also below. Note 
that {h{u)\n{u)) is purely imaginary (due to {n{u)\n{u)) = 1), so that this phase factor nullifies, if \n{u)) can be 
chosen to be real. An alternative representation of |(/'n(s)) is 

Ed 
, Ckn(s)\k(s))^ Cfcn(O) = 4n, (21) 

k—1 

Ckn{s) = 6kn + ecW (s) + e'cf^is) + e^clfl (s) + . . . (22) 



This fact also explains why in (14) it suffices to take a single initial wave-vector and not a superposition of them. Any superposition 
will bring in the adiabatic limit strong oscillations for non-diagonal elements of the resulting density matrix. This will reduce the 
superposition to the mixture of adiabatic cigen-vcctors, which amounts to studying the consequences of (14), and then tciking the 
average over the index n with certain time-independent weights. 
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Let us quote from Appendix B several basic formulas of the adiabatic perturbation theory: 



fc^n(s) = ^TT^t:^^ Afei(s) = Ek{s) - Eiis), (23) 



[1] _ (fc(s)|n(.s)) 
zA„fc(s) 

W= A^RinH")-(*K)l (25) 



where means that k = n is excluded from the summation X]fe=i- 

Eq. (23) for the first-order adiabatic correction is especially well-known [5, 6]. It is certainly less well-known that 

the consistent adiabatic perturbation theory generates another 0{e) term, i.e., ckU; see however [24]. This term is 
purely imaginary and it drops out from the lowest-order post-adiabatic corrections to the averaged force in (3). 

The representation (21) is gauge-covariant \(j)n{s)) e*""^*) |0„(s)) due to (23-25). Eq. (19) is also gauge-covariant 
IV'raC*)) ~^ e'""^°^|'^„(s)), and shows that the relative phase arg [(n(0)|^„(s))] is gauge-invariant. Hence it can be 
observed [25, 26]. Note that Berry phase factor in (19) is not gauge-invariant as such, apart of a certain specific case 
[25-27]. 

A. Compact expression for the non-adiabatic force 

Employing (17, 19) we now calculate 

F^ = {i>n\di,H\i,^) = (,/.„|a^jy|</.„) (26) 
= dMH\4>^ -'2^{d,MH\4>n) (27) 
= d^Er, + d^[{,t>n\H\ct>n) " E^] + 2e9(a^0„|</.„), (28) 

where 5ft and 3 mean, respectively, the real and imaginary parts. The factor S^E'n is the force generated by the 
adiabatic (Born-Oppenheimer) potential E^ = En{q). Thus the last two expressions in (28) represent the non- 
adiabatic force. We denote 

F, = i^oi + + e^i^^] ^ ^3^3] + _ _ _ ^ (29) 
where F^°' = d^En- In this context note from (21) that 

(<^„|F|.^„) -E^ = Afc„|cfc„|2 = 0{e^). (30) 

IV. FIRST- AND SECOND-ORDER POST-ADIABATIC FORCE. 

Using (22) and (28) we get for the first-order post-adiabatic force: 

fW = 29(5^n|n) = 2q^'^{d^n\do,n) , (31) 
where we always assume implicit summation from 1 to K over the repeated Greek (not Latin!) indices: 



A„Br 



= V (32) 

^ — ^a— 1 

Since ^{dij,n\dan) = — 9(9an|9^n), Eq. (31) leads to an efi'ective Lorentz [or gyroscopic] force [5]. 
Employing (21, 23) we obtain at the second order: 

= ^MEl^'^-l^nP + 23 {{d,n,\h) + {d,n\h,)) (33) 
= a,E;^^^|^ + 29(^(a,n|n,)+a,(n,|n)) . (34) 

Note that the non-local contribution (^nn drops out from (34), since both CnJ, and {n\dij,n) are purely imaginary. 
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Working out (33, 34) and taking the result together with (31), we obtain that in the first and second orders the 
averaged force can be generated by the following Lagrangian ^ 

-^[^'1(9,9) = yGa/3(g)9a9/3 + e^a(9)9«, (36) 
A„ = Q{dan\n), (37) 
= -2V 3?{(n|5^fc)(9<,fc|n)}. (38) 

Here is the vector potential that corresponds to the effective magnetic field (31) [5], and Ga/3 plays the role 
of a coordinate-dependent mass tensor, which is generated from (34) [6]. Note that Gap is a positive matrix, i.e., 
Gai34'a4>6 ^ 0; for any vector (f>a, provided that the the quantum system starts its evolution from the ground state: 
^fcn > 0. Now Gai3 Cannot be a positive matrix for all initial states of the quantum system, since, e.g., when d = 2 
(two level situation), one has Ga/j [excited state] = —Ga/j [ground state]. 

The complete classical Lagrangian to the second order is obtained by adding C^^'^^{q, q) and the Born-Oppenheimer 
potential En{q) to the initial (bare) classical Lagrangian Cq, given by (1): 



^2 = fEL(^) -V{q)-E^{q)+6^\q,q). (39) 



Note that when the time-scale separation is enforced by a large [bare] mass M of the classical particle, the post- 
adiabatic Lagrangian C^^'^\q,q) is small as compared to the large kinetic energy > m&ke this fact 

explicit, we rescale this kinetic energy to the slow time via e ~ 



A. Second-order post-adiabatic force, metric tensor and curvature. 

The kinetic part ^[M5afj + Gaii{q)]qaqp of the second-order Lagrangian (39) corresponds to a free particle moving 
on a Riemannian manifold with metric tensor [3]: 

gafiiq) = e^[MSc0 + G„/3(g)]. (40) 

There is an important particular case, where the complete Lagrangian (39) just reduces to this kinetic energy. This 
happens when i) the eigenvectors |n) can be chosen real, which then nullifies the vector potential Aa, ii) the bare 
potential V{q) and the Born-Oppcnheimcr potential En{q) compensate each other, V{q) + En{q) = 0, e.g., V{q) is 
zero from the outset, while En{q) turns to zero, since the eigenvalues of the quantum Hamiltonian H[q\ do not depend 
on the coordinates q (though the eigen- vectors do). 
Thus we focus on the purely kinetic Lagrangian 

\9c.p{qWe- (41) 

Once we are going to excersise on the Riemannian geometry, we recover for the velocities the explicitly contravariant 
notations [3] g". The metric tensor gaf} is then naturally covariant. The Lagrangian (41) yields the following equations 
of motion 

r + r^.g'^g^ = 0. (42) 



^ Let we are given a classical system with action /p ds£[g(s), ^(s) ], where £ is the Lagrangian, g is the vector of (generalized) coordinates, 

and 9 = gf • The Euler-Lagrange variational equations of motion ^ = 0, arc obtained when varying the action over the 

the coordinate-path g(s) assuming that the end-points are fixed; 6g(0) = 5q{S) = 0. This well-known set-up has a straightforward 
generalization for a Lagrangian C[q{s), q{s),q{s)] that depend on the acceleration [or more generally on higher-order derivatives of 

coordinates]. The corresponding Euler-Lagrange equations of motion read: ^ [^t] ~ ^' where now we assume that 

5g(0) = Sq{S) = Sq{0) = 6q{S) = 0. 
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This is the geodesic equation = 0, where the covariant differential of any vector C" is defined as 

DC"' = dC" + r«^C"^d<, (43) 
and where the connections F^^ relate to the metric tensor as [3]: 

= Ig"^ {9^9-^^^ + duQa^. - da9^.) . (44) 

Here g""^ is the inverse of the metic tensor: g°"^grT0 = S^, and where is the Kronecker delta-symbol. 

The first important question is wh(rt,her th(^ Riemannian manifold is curved or not. In the latter case it is possible 
to bring g^p to a diagonal and coordinate independent form by going to some new coordinates q' . The criterion of 
this is the Riemannian curvature tensor R'^^/j [3]. It will be more convenient to present the explicit formula for the 
covariant curvature tensor [3] 

-Ra/375 = 2 [dj^QaS + 5^55/37 " ^I^^CT " d^^Q^s] + g^" [^u0iTu,aS - ^i^0sTtJ.aj] , (45) 

where F^/j-y = jf^^F^^. Eq. (45) implies the following symmetry relations: 

Ra0iS = —RffajS = —RafiSj = R^Safi- (46) 

If the curvature tensor is non zero the manifold is curved. The manifold is not curved, if and only if R'^^^ = 0. Recall 
that for any vector C", the curvature tensor determines the non-commutativity degree of the covariant derivatives 
[3]: 



/~io. /~ia. /tar/a. /~ioi — ^ ( A1\ 



It is known that the curvature tensor determines the local behaviour of geodesies with respect to perturbing their 
initial conditions [3]. Let .7:"(.s, 0) be a family of geodesies, where .s is the time, and is a scalar continuous parameter 
which distinguishes the members of the family. Thus by the definition of the geodesic: 

Introduce a vector v"- = , which determines the deviation of two geodesies with slightly perturbed initial conditions. 
This vector satisfies the following Jacobi-Levi-Civita equation [3] ^: 

^ = R'^.^s^W. (49) 

The vector v" can be separated into two components v" = v^-^ + v^y. one orthogonal to u" {uaV^j = 0) and another 
one parallel to u". One can check with help of (46, 48) that the orthogonal component v^-^ satisfies the same equation 

(49), while the parallel component satisfies the geodesic equation (48). 

Below we calculate the curvature for the simplest example of two classical coordinates and g^. The fact of having 
only two coordinates simplifies the formulas for the curvature. Eqs. (46) imply that there is only one independent 

component of the [covariant] curvature tensor, which can be chosen to be i?i2i2- AH other components are either zero 
or equal to ±i?i2i2. One can check that now Rap-yS is expressed as 

R 

Ra0'yS ^ -^[ga^gpS - 9a59l3j], (50) 

D <■>-, f:IS r, 2iii2l2 , . 

R = 9 '9 Ral3-/S = 2~' (51) 

511522 - 512 



* To derive Eq. (49) note that the very definitions of and v°' imply dpu" = u^dpv°', which amounts to u"^?;'' = v^^u^^. Now 
calculate directly ^, ^ recalling (47) and noting that u^gU^ = due to (48). 
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where R is the scalar curvature. The latter thus determines the whole curvature tensor for the present two-dimensional 
situation. Substituting (50) into (49) and recalling that one can take UaV" = in this equation, we get 



(52) 



Note that Ujju'^ does not depend on s; see (48). 



We now set to calculate the curvature tensor -RJ^q,^ for the simplest possible example, where there are only two 

classical coordinates , and the quantum system has only two energy levels. For further simplicity we assume 
that the quantum Hamiltonian is real. This means that the Hamiltonian is a linear combination of the first and third 
Pauli matrices: 



H 



(53) 



The eigenvalues and eigenvectors of 11 read respectively 



E_ 
1+) 

I-) 



1 



P-T 



(54) 
(55) 

(56) 
(57) 



It is seen that the adiabatic energies Ej^ and E- cross at p = 0. 

We shall study in separate the case when the quantum system starts at t = from its ground state |— ), and from 
the excited state |+). 



1. The ground state. 



The metric reads form (40) and (54-57): 



511 



M + 



4p5 



gi2 = 921 = -e 



4p5 



The determinant and trace of the metric read 

det[ff] = e^M ( M+ ^ 



1 



4p3 



522 



M + 



4p5 



tr[5] 



1 



(58) 

(59) 
(60) 



(61) 



It is seen from (58 61) that both the determinant and the trace of ga3 are positive; thus the eigenvalues are positive 
as well. This situation refers to a usual classical mechanical particle, which is enforced to move on a two-dimensional 
surface. For the scalar curvature we get from (44, 45, 51) and (58 61) 



R 



3(1 + 16Mp3) 
2e2Mp2 (1 + AMp^f 



(62) 



R is strictly negative. Returning to (52) we see that since the metric is positively defined [see (58-61)] u°'Ua is 
always non-negative. Then the negativity of R in (62) implies that the geodesies are unstable with respect to small 
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perturbation of initial conditions, because (52) corresponds to a harmonic oscillator with an inverted (though space- 
dependent) frequency This unstability might have implications for the validity of the adopted adiabatic assumption, 
a question which we plan to study elsewhere. Note as well that R is singular at p = 0, where the adiabatic energy 
levels cross. 



2. The excited initial state. 



Now we assume that the two-level quantmn system starts its evolution from the excited state |+). This case leads 
to more interesting possibilities, since now the metric reads: 



511 = e 



M- 



512 = 521 = e 



4p5 
4p5 ' 



522 = e 



M- 



4p5 



Hence the determinate and trace of g read, respectively. 



det[g] 



M- 



1 



tr[5] 



M 



4p3 



Since the metric (63, 65) relates to (58, 60) with transformation M 
scalar curvature directly from (62) 



-M and e ^ ie («^ 



3(16Mp3 - 1) 
2e2Mp2 (1 _ 4Mp3)2 



(63) 

(64) 
(65) 

(66) 

-1), we get for the 
(67) 



When the particle moves sufficiently far from the origin = q'^ = (i.e., when 4Mp^ > 1), the metric is positively 
defined and the curvature is positive. According to (52) this means that the geodesies are not sensitive to perturbations 
in initial conditions. At AMp^ = 1 the metric tensor changes its signature, so that for 4:Mp^ < 1 it has one positive 
and one negative eigenvalue. At 4Mp^ = 1 the scalar curvature is singular. We expect that the adiabatic assumption 
will become problematic in the vicinity of the singularity, but the possibility for the particle to "tunnel" between 
subspaces of different signature is striking and deserves a more detailed investigation. 

Since the metric tensor is not positively defined for 4Mp^ < 1, (52, 67) show that for | < 4Mp^ < 1 the geodesies 
with initial condition Uau" < can become unstable 

For even smaller values WMp"^ < 1 of p the curvature becomes negative. Now the unstable geodesies have UaU" > 0, 
while those with Uau" < are (at least locally) stable. 

It is thus seen that the initial ground versus the excited state of the quantum system produce rather different 
dynamic behaviour for the classical system. 



V. THIRD-ORDER POST-ADIABATIC FORCE. 



We now turn to studying the post-adiabatic force at the order e^. The calculations here are more involved, 
though their general pattern employing the adiabatic perturbation theory and then reconstructing the effective 
Lagrangian — remains the same. The calculation details being presented in Appendix C, we shall quote the main 
result: At the order the non-adiabatic force acting on the classical system is still Lagrangian [see Footnote 7]: 

d d _ d^_a d 

^ ds dq^ ds^ dq^ dq^ 

C.^%q,qA] = [faff-yQaq/SQi - Zal3(iaQ0] , (69) 



^^i^^' = (68) 



Such a local instability leads to chaos, if the {q^ , g^)-manifold is compact. This is not the case for the considered situation, though it is 
presumably not very difficult to compactify the manifold, keeping the conclusion on the local instability of geodesies. 
In the General Theory of Relativity Uau" < is prohibited by causality; for massive particles Uau" > and can be normalized to 
Uau" = 1, while for photons Uau" =0 [3]. Ifowcvcr, for the present classical theory with a well-defined global time s nothing prohibits 
us to consider the class of geodesies with UaU°' < 0. 
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where stands for the third-order Lagrangian, and where we defined 

/a/37 = Q[{n\d^n){Nf,\N„) + {d^Ni,\N^)], (70) 
ZaIS = ^[{NplN^)], (71) 

It is seen that besides the expected third-order polynomial over the velocities fafS-yQaQfjQ'y, the third-order Lagrangian 
£P1 contains a linear dependence on the accelerations ija- The corresponding coupling matrix Za/siq) has to be 
antisymmetric Za/siq) = —zpaiq), since any term (jya/sQaqp with a symmetric (j)c,f3 = (p/Ba, can be reduced (up to a total 
differential in time) to a third-order polynomial over the velocities. 

The total Lagrangian describing the classical system including the three-order terms will include the previous order 
non-adiabatic forces and the bare classical Lagrangian, Cslq, q, q] = (39) + £[^1 [q, q, q], or 



e 



^sb, q-, q\ -y{q) - En{q) + e^a (<?)<?« + -^[MS^fj + GalsiqMocqp + e [/a/379a3/397 - Zaisqaqi}] , (73) 

while the equations of motion it generates is [see Footnote 7] 

d d d d 



^dsdqij, ds^ dq^ dq^J ^^^'^'^'^^ ^' ''^^^ 
These equations of motion contain third-order time-derivatives qa^ of coordinates, i.e., they can be written as 

2e^Zafjqf^ = (Paiq, q, q) . (75) 

Thus when det[za/3] 7^ — and this is gcnerically the case for even number of classical coordinates — the third- 
derivatives can be expressed through {q, q, q). This means that the dynamics described by (74) needs three independent 
(sets of) initial conditions at the initial (slow) time Si: 

(g(,si), g(si), g(si)). (76) 

For an odd number K of classical coordinates, one has det[2;a/3] = 0, since Zaff is anti-symmetric. Generically, the 

matrix Zap will have only one eigenvalue equal to zero. Let us denote this eigenvector by ya\ Va^Zap — 0, and let 

Ua^ with 7 = 1,. ..,K — 1 be the eigenvector of Zap with non-zero eigenvalues X^'^h Eq. (75) produces 

2e3AW y^^\f = y^jWiq, q, q), for 7 = ■ ■ ■ , K - 1, (77) 

= ?/L°Va(9,g,9)- (78) 

Now the initial conditions {q{si), q{si), q{si) ) at the initial time Si cannot be anymore taken independently from each 
other, because (78) imposes a constraint on them. Provided that {q{si), (j(si), q{si) ) satisfy this constraint, (77) gives 

K — 1 equations for components of ga • Another equation for components of qa can be obtained by differentiating 
(78) over time t and taking t ^ U. 

Before closing this discussion on the initial conditions let us note the following aspect. The quantum-classical 
equations (2, 3) have the following well-defined initial conditions at the initial moment t = of the fast time t: 
1^1/(0)), g(0) and g(0). On the other hand, as we saw above, the autonomous classical dynamics starts to depend on 
higher- derivatives of the coordinate(s). The reason of this diflFerence is that the initial conditions of the autonomous 
classical dynamics are to be imposed at an initial value Sj = et-^ of the slow time, where t\ should be still sizably 
larger than t = Q. The difference between the original initial conditions of the slow variables and their effective initial 
conditions after eliminating the fast variables is known as the initial slip problem. It is well recognized in theories 
dealing with elimination of fast variables [35] . There also exist more or less regular procedures of relating the original 
initial conditions to effective ones [35]. We shall not dwell into this issue anymore, because in the present paper we 
are interested by autonomous classical dynamics for sufficiently large (fast) times, where the precise relation with the 
original initial conditions is not directly relevant. 



The construction described around (77, 78) is conceptually not very diflfcrcnt from its simplest analog: Consider two classical degrees of 

freedom with coordinates x and q. Let the corresponding Lagrangian be ^ V{q^x). Note that this Lagrangian does not contain the 

kinetic energy for the x-particle, i.e., the kinetic energy matrix is degenerate. The Lagrange equations of motion read: q = —Vq(q,x) 
and V^{q, x) = 0. The second equation is a constraint on admissible values of q and x at any time. In particular, it confines their initial 
values. Now the initial conditions amount to q{si), g(si) and x(si) provided that the constraint is satisfied. One is not free in choosing 
the initial velocity x{si). The latter is determined from differentiating the constraint over time s and taking s — > Si. 
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A. Kinematics. 



The dependence of Cz[q^ q, q] on accelerations implies conceptual changes in the kinematics of the classical system, 

as we now proceed to show. 

First we note that the momentum of the classical system is defined via the response of £3 to an infinitesimal 
coordinate shift — > + Sq^, where Sq^ does not depend on time [3]: 



dCs d 
dq/j, "^As 



5£3 d 5£3 



dq^ ds dq^ 



(79) 



where we used (74). Thus the momentum is defined as 

9/:3 d 5£3 



dqf. As dq^ 



(80) 



implying that the equations of motion can be written as = 

If £3 would not depend on (which is generically not the case), the corresponding momentum will be conserved 
in time. Note that consists of the usual part and the anomalous part — ^fy^ that comes solely from the 
dependence of the Lagrangian on the acceleration. Using (73) we get for the momentum 



(81) 



where /^'J™'' = ^Sn/n[a/37] is the completely symmetrized expression (70); the sum is taken over all six permutations 
n of three elements. It is seen that the expression for the momentum does depend linearly on the acceleration. One 
half of the acceleration-dependence comes from usual part , while another half comes through the anomalous part 

resulting altogether in 2e^z^aqa in (81). 
The energy corresponding to the Lagrangian C3[q,q,q] is obtained via looking at the total time-derivative of 
Jl^3[q,q,q\- 



-T:^3[q,q,q\ = ^ — q^ + -^Qfi + -^-r-^' 

dt dq^ dq^ dq^ ds^ 



(82) 



where we noted that C^lq, q, q\ does not have any explicit time-dependence. Employing equations of motion p^ = 
(82) results into energy conservation: 



d^ 



E = p^q/j. + -Q^% - l-'-i- 



For our case the energy E reads 



E = — [MSa/3 + Gaffjqaqff + '^e^ fal3iqaq/3qj + 2e^ + V{q)+Er,{q). 



(83) 



(84) 



Note that the vector-potential Aa{q) cxpectedly drops out from the expression of energy [3]. However, the acceleration- 
dependent part of the Lagrangian does contribute directly into the energy. In fact, the whole third-order Lagrangian 
(69) is multiplied by 2 and enters into the energy. 

Let us now turn to the angular momentum tensor, which is defined via the response of £3 to an infinitesimal 
rotation (i.e., a distance conserving linear transformation) [3]: q^ +bj^„5qrj, where oj^^ = —(jJo-h- 



SC. 



dq, 



-qp 



dC3 

dqa 



qp, 



Paqp 



-^qp 

dqa 



(85) 



where we again used (74). The full momentum tensor is now defined as [recalling w^^- = —LOati]'- 

Map = Paqp - Ppqa + -q^1i3 " '^^^ 



ap ■ 



^ap, 



(86) 

(87) 
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so that when £3 is rotationally invariant, Map is conserved in time. One part of this tensor is the usual orbital 
momentum L^ij = Palp ~ PpQa- The remainder — non-orbital momentum, or spin — arises due to the dependence of 
the Lagrangian on the accelerations, and it is a second-order polynomial over the velocities: 

= e^[z/3jq^qa - ZajQ^qp]. (89) 

In the simplest two-coordinate situation S12 = —S21 = £^-221(91 +91)) which means that the spin tensor is proportional 
to the velocity square. 



1. Zitterbewegung effect. 
Now note from (80, 88, 89) and from zpa = —Zap that the momentum can be written as 



- dq. + dt 



r = mf3, (90) 



which means that the anomalous part P/j, — of the momentum is driven by the time-derivative of the velocity- 
projected spin-tensor. 

An expression similar to (90) — relating the momentum to the projected time-derivative of the spin — appears in the 
(relativistic) Dirac electron theory; see [14] for a review. There the fact that the total angular momentum is a sum 
of the orbital part and the spin part, as well as the fact that the velocity and the momcmcntum arc different objects 
and are not simply proportional to each other via the mass, arc the consequence of the relativistic invariance for the 
electron. The very effect of the spin time-derivative contributing into the momentum was named zitterbewegung, since 
for the free Dirac electron this contribution brings in an additional oscilatory motion [14] . In a more recent literature, 
the zitterbewegung effect is also derived via Lagrangians contaning the higher-order derivatives of coordinates [28, 29]. 

There arc, however, several aspects that distinguish (90) from the zitterbewegung effects already known in literature. 
First, we do not have a relativistic invariant theory; for us relation (90) emerges due to the fact that the classical 
system is open. Second, we do not have to have the conservation of momentum and of angular momentum for deriving 
(90). Both these quantities are generically non-conserved in our situation (ultimately since the system is open), but 
relation (90) still holds generally due to the specific, anti-symmetric form (69) of the acceleration-dependent part of 
the Lagrangian. 

We close this part by re- iterating its main findings: due to interaction with the fast quantum system the classical 
system gets a spin [non-orbital angular momentum], which is related to its momentum via zitterbewegung effect. 



B. Hamiltonian description. 



Further insight into the structure of the effective classical dynamics is gained by studying its Hamiltonian description. 
Within the order the Hamiltonian description is straightforward. However, the third-order dynamics has a non- 
trivial Hamiltonian structure, as seen below. Rewrite (73) as 

■^3(9, q, q) = L3{q, q) - e^ZafsqaQp, (91) 

where the higher- derivative term is explicitly separated out. Instead of (91) we now introduce the following extended 
Lagrangian: 

L[q, V, n] = Waiqa - Va) - e^ZapVaVfS + Lsiq, v), (92) 

which is a function of three set of variables: q = (gi, . . . , qx), v = (wi, . . . , vk) and tt = (tti, . . . , ttk)- It should be 
clear that if we treat q, v and tt as coordinates, then the Lagrange equations generated by L[q,v,n] are equivalent 
to those generated by Cz{q,q,q). At this point tt is considered as a part of the overall set of coordinates. It may be 
equivalently viewed as Lagrange multipliers. If L[q, v, n] were not depend on v — that is £3(9, q, q) were not depend on 
ij- we would write the velocities v = v{q, n) as functions of the coordinates and momenta, and end up with the usual 
Hamiltonian description with q and tt being the canonical coordinates and momenta, respectively. Though L[q,v,iT] 
does depend on w, it can be still Hamiltonized following to the method advocated in [31]. 
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Once the triple q, v, w is considered as coordinates, we introduce a separate notation for it 

Q = {Ql,...,Q3K) = {qa,Va,na). (93) 

Now the Lagrangian (92) reads 

L[Q] = Aa[Q]Qa + n[Q], n = L3[Q], (94) 

where the index a runs from 1 to 3K, and where Aa = (tTc, e'^2^at'/3) 0) is deduced from (92). As we show below, Ti-IQ] 
will play the role of Hamiltonian. Eq. (94) generates the following Lagrange equations of motion: 

^ab{Q)Qb = ^, (95) 

In block-matrix notations f2 reads 

/ Y I\ 

(97) 

y 

where each element in the above matrix \s K y. K matrix: 

Yap = ^^V-ydaZjff, Za/S = -2e^ Za0, laj} = (98) 

and where / is the K x K unit matrix. Provided Z is invertible, the inverse of O reads [block-matrix notations] 



fi-i = - Z-iy^ . (99) 






Y 


I 


-Y^ 


Z 





-I 










For an even K the matrix Z is gcncrically invertible; compare with our discussion after (74). Since Qab is invertible, 
antisymmetric, closed and it satisfies (95), Qab defines a symplectic structure [22]. Then Tl[Q] plays the role of 
Hamiltonian. These two ingredients are necessary and sufficient for the Hamiltonian description [22]. In particular, 
for any two functions C{Q) and D{Q) the Poisson bracket is defined as 

f)(^ f)r) 

{CiQ),DmPB = n-:, — — . (100) 
The equations of motion (95) are now written as 

Qa = {Qa,n[Q]}pB. (101) 

Note that the Poisson brackets are non-linear. It is seen from (97, 100) that Za0 and its derivatives define a non-trivial 
symplectic structure for the system. 

The matrix Z is not invertible for an odd K. The Hamiltonian description in this case is still possible, but requires 
more care in explicitly accounting for constraints; compare with our discussion after (74). 

VI. THE FOURTH ORDER LAGRANGIAN. 

Here we briefly report on the fourth-order Lagrangian. Since the calculations now become very complicated, we shall 
restrict ourselves to the situation where the classical system has just one single coordiantc q. For further simplicity we 
assume the quantum system has real adiabatic eigenstates. In fact, the main purpose of this section is to illustrate that 
the fourth-order Lagrangian again depends linearly on the highest-order time-derivatives of the classical coordinate. 



1^ Closed means that ^^f2„i, + j^^Qca + g^^bc = 0- This feature is automatic from the definition (96), and it ensures that the Poisson 
brackets defined via Uab does not change in time [22] . 
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Following the same lines of calculation for the thrid order non-adiabatic force presented in AppendixC, and assuming 
a single coordinate classical system and real eigenstates for the quantum system, the non-adiabatic force acting on 
the classical system in the fourth order is described by the following Lagrangian 

\ds^dq(^) ds^dq^dsdq Oq J ^^'^'^'^ ^' ^ ' 

[q, q, q, 9^'^] = [og^ + bqq'' + wqq^'''^] , (103) 

where q^^^ stands for the third order time derivative of q, and where >CW represents the fourth-order Lagrangian. Note 
that the dependence on the higher-order time derivatives q and q^^^ is linear. The coefBcients a, b, and d are given as 

nk 

where \N) is given by (72): \N) = J2[ih^\k}. 

Then the total Lagrangian describing tSe one dimensional classical system reads 

q, q, q^'H = -V{q) - i?„(g) + ^(M -h G)q' + (aq^ + bqq' + wqq^^^) , (107) 

where the first and the third order terms vanish due to the assumption of real eigenstates of the quantum system, 
and where a, b, and w are given by (104)- (106) and G is defined as 

G=-2y ^ 'f ' ' . (108) 
The kinematics of this Lagrangian can be developed along the same lines as in the previous section. 



VII. SUMMARY 



We have studied the post-adiabatic equations of motion for a slow classical system which is coupled to a fast 
quantum system. The slow versus fast motion is controlled by a small ratio e of the characteristic times. The general 
problem we addressed is to find an efi^ective Lagrangian that describes the dynamics of the classical system. The 
following facts were known for this problem: (1.) In the order the effective Lagrangian differs from the bare one 
by the; Born-Oppcnhcimer potential energy. (2.) Berry and Robbins have shown that the e-force term in the effective 
Lagrangian corresponds to a magnetic field, which is related to the geometric phase [5]. (3.) Weigert and Littlejohn 
and Goldhaber has recently shown that in the order the effective Lagrangian gets an additional kinetic energy term, 
which is a second-order polynomial over the classical velocities [6, 7]. 

In this work we obtain the following new results. 

i) The post-adiabatic reaction force is proved to be Lagrangian up to the fifth order in e. We conjecture that at 
every order of e the effective dynamics of the classical system can be derived form a Lagrangian. 

a) Within the order e"^ the effective Lagrangian linearly depends on the accelerations of the classical system. 

We argued that this result is important, because it provides a physically well-motivated scenario for the emergence 
of higher-derivative Lagrangians for open classical systems. This result should be contrasted to the usual open-system 
approaches, which can also produce forces depending on higher-order derivatives (e.g., the Abraham-Lorentz force 
in electrodynamics), but those forces are dissipative (non-Lagrangian). We also explained that this result does not 
depend on the quantum character of the fast system, and will be present as well if the fast system is classical and 
integrable [in the sense of [22]]. 

The presence of higher-derivative terms can be tested by essential influences they bring on the kinematics of the 
system. First, they modify initial conditions; in our case this means that the trajectory of the classical system on the 
slow (coarse-grained) time starts to depend on acceleration; see our discussion around (75-78). Second, the conserved 
energy of the slow classical motion does depend on the acceleration. And third, the presence of higher-derivative 
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terms naturally separates the total angular momentum into the sum of orbital momentum and spin. We show that 
this spin satisfies an exact analogue of the zitterbewegung relation. 

in) We also found interesting results for the classical autonomous dynamics within the order e^, where the motion 
generated by the effective classical Lagrangian can be mapped to geodesic curves on a suitable Riemannian manifold. 
Operating with the simplest possible example — two classical coordinates interacting with a two-level quantum sys- 
tem — we show that the Riemannian manifold is essentially curved solely due to the kinetic energy generated by the 
fast quantum system. The scalar curvature is frequently negative indicating that the classical trajectories [geodesic 
curves] are unstable with respect to small variations of the initial conditions. The metric tensor generated by the fast 
quantum system can change its signature as a function of the two coordinates. Physically this means a transition 
from an Euclidean to pseudo-Euclidean manifold, emergence of a time-like coordinate and etc. This result deserves a 
careful elaboration. 
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APPENDIX A: ON THE DERIVATION OF THE QUANTUM-CLASSICAL DYNAMICS 

The argument follows basically to [19, 32]. Let we are given a two-degrees-of-freedom quantum system with 

2 2 

Hamiltonian ^ + V{q) + H{q,x) + j^, where q and x are operator coordinates, and where p and tt are operator 
momenta. For simphcity we shall assume that the initial state is factorized over these two degrees of freedom. 
The Heisenberg equation generated by this Hamiltonian read: 

^-P- (Al) 

f^=-V'{q)-H'^{q,x), (A2) 

I = (A3) 
at m 



-K{q,x). (A4) 



d-TT 

Now the motion of the (p, q) degree of freedom is separated into two part: 

p{t) = p{t) + pf, q{t) = q{t) + qf, (A5) 

where p{t) and q{t) are the averages over the initial state, and where pf and Xf are the fluctuations. By definition 
these operators satisfy 

Pf = qf = 0. (A6) 
We now substitute (A5) into (Al, A2, A4) and expand (A2, A4) over the small qf. 



dq dqf _ p pf 
dt dt M M' 



(A7) 



^ = -HUq, x) - H':^{q, x) qf + 0{q}). (A9) 
Averaging these equations over the initial state we obtain 



(AlO) 



dg ^ p_ 
dt M' 

g = -VW) - - H'^M ^)lf + (All) 

Htt 

— = -Hl^{q,x) - H'^^{q,x)qf+0{q)). (A12) 



18 



If in (All, A12) the terms proportional to 0{qf) arc neglected we get into a quantum-classical equations, where (p.q) 
is the classical degree of freedom. If, however, there are physical reasons to expect that the initial state will remain 
factorized over the considered range of times [sec [10] for a detailed discussion of such situations], then we have to 
neglect only terms 0{q'j), because the terms 0{qf) drop out due to (A6). 

APPENDIX B: ADIABATIC PERTURBATION THEORY. 

Here we outline the adiabatic perturbation theory as developed in [23, 24]. This appendix is completely self- 
contained and can be read independently from the main text. 
Consider the non-stationary Schrodinger equation 

zdt\^it))=H{et)\^{t)), (Bl) 

where = 1, e <C 1 is a small parameter, and where H{et) is the slowly changing Hamiltonian. Let us define 

s = et, (B2) 

for the slow time. Let En{s) and |n(s)) be the adiabatic eigen-energies and eigenvectors, respectively: 

H{s)\n{s)) = Er,{s)\n{s)), {n{s)\m{s)) = Snm, n=l,...,d. (B3) 

We assume that the adiabatic energy spectrum is non-degenerate for all s. Let us also assume for simplicity that the 
system starts from the initial state equal to one of the adiabatic eigen-states: 

|*(0)) = |n(0)). (B4) 

The first step in any adiabatic approach is to separate the slowly changing quantities from the fast ones. To this 
end we look for the solution of (Bl) as 

|*(e,f)) = |^„(e,s))e*""W, a^it) = - f dt Ei{ei), (B5) 



JO 

where an{t) is the dynamical phase. It is clear that e'""^*) changes fast, i.e., as ~ e'/^, for slow times s. Putting (B5) 
into (Bl) we get 

ie\Me, s)) = [H{s) - £„(s)] |Vn(e, s)), (B6) 
where dot means differentiation over s. Now we expand the slow wave-function, Vn(e, s), over e 

|^„(e, s)) = e^T"W [ \n{s)} + e\n,{s)) + e'\n2{s)) + ...] , (B7) 

where 

jn{s)=i f du{n{s)\n{u)), (B8) 



is the Berry phase factor. We separated it out to facilitate further calculations. 

Substituting power series expansion (B7) into (B6) and comparing terms of equal order of e, we get 

0={H-En)\n), (B9) 
i\h) - i{n\n) \n) = {H - E^) jni), (BIO) 
ijni) - i{n\h) Im) = {H- En) (na), (Bll) 



Eq. (B9) holds automatically. To solve the higher order equations we introduce the projection operators P and Q: 

P=\n){n\, = ^Jfc)(fc|, (B12) 
P + Q = l, PQ = QP = 0. (B13) 
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where means the term k = n is excluded from the summation X^f—i- 
Operating Q from left to the both sides of (BIO), we get 

i^Jfc)(fc|n) = Y,^Akn\k){k\ni), Afe„ = Ek - En- (B14) 

Since Afc^„ is non-zero [energy levels are not degenerate], we get 



(fc|nx) ^ cW „ = -if^. (B15) 



This determines Q\n\), but we still have to find P\n\). Let us define 

\n,) = c\^i\n) + \ni), (B16) 

cW ^ (n|ni), (B18) 

where c[^i has to be found. To this end we multiply both sides of (Bll) from left by P [recall that (Bll) is obtained 
from term in expansion (B7)] 



where jni) is found from (B16): 



Using (B19) we get 



(njni) = (njni) (njn), (B19) 
\ni)=c^l\n) + c^^i\h) + \hi). (B20) 

c^l{s) = -{n{s)\ni{s)), (B21) 



which together with (B4) solves as 

cW (.) = - f Mn{u)\ni{u)) = -.V ' f du '^'^"^'"/f (B22) 

Jo Jo 'Ank[U) 

It is seen that c[m(,s) is purely imaginary. 

Analogous argument gives the higher order non-adiabatic corrections \nm{s)), (m > 1) in (B7) 



^)=E/4'::'i^)+ct^w, (B23) 

n 



[m] _ ci"„^' - i{k\hm-i) 



= -E//'dncI'::](«) {n{u)\k{u)). (B25) 

Altogether |?/'„(e,s)) in (B7) can be written as 

|V„(e,s))=e^^''(^)^^Cfe„|A;), (B26) 

Some relations between coefficients c/j„ can be uncovered without knowing their explicit form, but rather looking 
at the normalization condition: 

Ejcnfcl' = 1, (B28) 
which should be satisfied at each order of e. For the first two orders 

|c„„|2 = 1 + 2e5ft{cW } + 2e^5ft{c[2] } + e^jc™ + ^(e^), (B29) 
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which brings in the following two relations in the orders e and e^, respectively, 

25R{cPl } + |cW |2 + Y^y^i? = 25?{ca} + (mlm) = 0. 

Let us work out explicitly the second-order coefficient c|.^^^(s). Using in (B27), {k\h\) = ■j^[{k\ni) ] — {k\ni) we get 



(B30) 
(B31) 



cgj.) = cW j5)cW(a) + ^-j^ [cW^ ((n|n)-(fc|fc)) [cW„] + ^^^^^ 



(B32) 



1. Precision of the adiabatic approximation 

An important problem of precision of the adiabatic approximation was studied in [23, 33, 34]. Simplifying previous 
results on this problem, Hagedorn and Joye have proven the following result [23]. Following to (B7) let us define 



|V^(e,s))=e'^"(^)^e'=|nfe(s)). 



(B33) 



fc=0 



Let {a} define the integer part of a real number a, and let we are given a positive number g. Then it is shown that 
[23]: 



\^{ah}[e, s)) - |Vexact(e, s)) < C(ff)e-r(«)/^ 



(B34) 



where jV'cxactle, s)) is the exact solution of the original equation (B6) ['(/''f^/'^ (e, 0)) = |'(/'exact(ej 0)) at the initial time, 
and where C{g) and T{g) are bounded positive functions of g. This results implies that the precision of the adiabatic 
approximation is exponential over e. 

APPENDIX C: CALCULATION OF THE THIRD-ORDER POST-ADIABATIC FORCE. 

Within the present Appendix means X]fc=i k^ni ^^'^ A = A' = where s is the slow time. 
The third-order post-adiabatic force is given by (28), where for — En we should employ (30) and then 

(22), while for 3(9^0„|(/)Jj) we directly use (20). Having done these, we select terms a and end up with 



~2~ 



ds 



(CI) 

(C2) 



[21 f [2] 

where wc additionally employed (B15) and |n2) = c„„|n) + X^/jCj.^|fc) when going from (CI) to (C2). 
Looking at (B32) we introduce the following notation 

.[2] 



_ „[!] „[11 , ^[2] 



Using (B31) we exclude KcJ^it and get 



~2~ 



-'^[{d^n\n')] (nilm) - ^9 [(5^n|n)] ^(ni|ni) - [(n|n')] d^{ni\ni) 



(C3) 

(C4) 
(C5) 



Let us denote: 



ni)=c^l\n) + \ni). 



(C6) 
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Using this notation and expanding 3(9^ni|n'j^) wc get after some algebraic steps that all non-local terms in (C4, C5) 
are cancelled out. This means that (C4, C5) are written as 



pl3] 1 A 1 



Using (B32) (C7, C8) can be represented as 



ds 



" ^nk 



After some algebra one can show that (C9) can be generated by a Lagrangian 



where 



(C7) 

(C8) 

(C9) 
(CIO) 

(Cll) 



Wt = l^{{n\d^ri)) mNaW^)] + ^9((n|5„n)) 5i[(7V^|iV^)] + l^{{n\d^n)) ^[{N^\N^)]. (C12) 



Here we defined 



^nk 



(C13) 



Recall that an implicit summation is carried over indices a, f3 and 7. Note that ha^-y is symmetric with respect to 
any permutation of indices a, (3 and 7, so that 



We now work out (CIO): Then 



y' (-— 

^nk 



> / 1 



"fen 



^E;4'(^'IOj(S.n|fc) 
1 ^'y\iuv 



■3 [(Af^|A^a)9a + 9ag/3(iVM|S/3Afa)] 



(C14) 

(C15) 
(C16) 
(C17) 



Combining together (C17, C18) we get 



(CIS) 



(CIO) = gag/39(5^iV/3|iV«) + g„g^q^9(a^iV;3|a7iV«) - —9 {{N^\N^)q^ + (?ag/3(iV^|5/3iV„)] (C19) 

= 9(iVa|7V^) + ga?/39 [a;3(Ar,|iV^) + (a^jJV^IAT^) + (a„7V/3|Ar^) + (a^A^/sliVa) ] (C20) 

+ gag0g7$>[a^(a;37Va|Ar^) + (9^Af/3|5^iVa)] (C21) 

We define 



(C22) 
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and write 



(C20) + (C21) 



d2__9_ _ d__d_ ^ _d_ 

ds^ dq^ ds dq^ 



Zapqaqfj 



(C23) 



d^{d0N„\N^) + {di,Np\d^N^) - -dlJNp\N^) 



where in obtaining (C24) we employed: 



(C25) 

(C26) 
(C27) 
(C28) 



The quantity inside of the square brackets in (C24) is symmetric with respect to permutation of indices /i, a and /3. 
We now define 

= [ {d0N^\N„) + (d^N^lN^) + (d^N^lNp) + {dc,N0\N^) + {d^N^lNfs) + {d^N0\N^) ] (C29) 



and obtain 



(C24) + (C25) 



d__d d_ 

dsdq^ dq^ 



(C30) 

(C31) 
(C32) 
(C33) 



Working out (C32, C33) we finally obtain: 

(C32) + (C33) = 0. (C34) 
Thus the third-order post-adiabatic force is purely Lagrangian (though containing higher-order derivatives): 



~2~ 



da d^ d d 
ds dq^ ds^ dq^j, dq^ 



(C35) 



